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A model-dependent analysis of the hadronic substructure in t~ — > 7r~7r~7r + i/ T decay is reported. The decay is 
dominated by the process r ~~ — » a^(Yl&S)v T , in which a\ (1260) — > p°7r~ in an S-wave decay. Amplitudes involving 
ai(1260) decays into isoscalars, especially ai(1260)~ — > /o(600)7r~, are large. t~ — > 7r~7r~7r + ;y T decays via the 
pseudoscalar 7r~(1300) are small. These results support the resonant substructure reported in the previously 
reported analysis of r~ — > ty°ty°tt^u t decay mode from the same CLEO II data sample. 
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1. INTRODUCTION 

The t~ — > (37r) _ f r decay has been the sub- 
ject of much interest over the years. Because of 
the transformation properties of the weak current 
under parity and G-parity, r-lepton decay to an 
odd number of pions is expected to occur exclu- 
sively through the axial vector current, ignoring 
isospin- violating effects. Since J p , and I G are 
good quantum numbers of the hadronic current 
from weak decays, the weak axial current pro- 
duces systems of pions with J p = 0~ or 1 + , odd 
G-parity and odd numbers of pions. 

The (37r) _ system in t~ — > (3ir)~is T is ex- 
pected to be produced dominantly through the 
poorly understood J p = 1 + a^(1260) where 
0^(1260) — > p°n~ in an S- wave decay. Although 
it is possible that the (3n)~ system can be pro- 
duced through the J p = 0~ 7r _ (1300), this is 
expected to be suppressed by the PCAC (par- 
tially conserved axial current) hypothesis. It has 
been known for some time that simple models 
do not provide a satisfactory description. For 
example, ARGUS W found that the model of 
Kiihn and Santamaria [ 2 ] and that of Isgur, Morn- 
ingstar, and Reader [3 did not describe their 
t~ — > (37r)~ v t data satisfactorily. 

Here we report the preliminary results of a 
model-dependent fit of the substructure in r~ — > 
7r~7r~7r + ^ T decays. The work is based on the 



analysis of ~ 145,000 r 



7r 7r t: + i>t events 



"representing the CLEO Collaboration. 



0. In addition to fitting the substructure with 
various a^(1260) decay modes, fits for the pro- 
duction of the (37r)~ system via the J p = 0~ 
7r~(1300) are included. These results are com- 
pared to the substructure in t~ — > 7r 7r°Tr _ ^ r 
previously reported ||. 

Throughout this paper the charge conjugate 
versions of decays are implied for brevity and clar- 
ity. 

2. EVENT SELECTION 

This analysis is based on 4.67 fb _1 of e + e~ col- 
liding beam data collected at center-of-mass en- 
ergy energies E cm ~ 10.6 GeV with the CLEO II 
detector at the Cornell Electron Storage Ring 
(CESR). This corresponds to 4.3 x 10 6 e+e" -> 
t + t" events, — > Tr~7r~ir + is T events are se- 
lected from a topology in which the r + decays 
into a single charged particle plus a v T and the 
r~ decays into three charged particles plus a v T 
(1 versus 3 charged-particle topology). Selected 
events are required to have four well-constructed 
charged tracks coming from the interaction re- 
gion, and the net charge of the tracks is required 
to be zero. No more than one track is allowed to 
be identified as an e~ or e + . The most isolated 
track is identified, and the remaining three tracks 
are required to be > 90° away from it. Events 
with photons (isolated electromagnetic showers 
not associated with a track) are rejected. Events 
that satisfy the hypothesis 0.485 < M T + V - < 
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0.510 GeV are rejected in order to eliminate 
events containing — > tt + tt~ decays. Finally, 
the 2-dimensional cut shown in Figure |l| is used 
to reject events with missing tracks and/or en- 
ergy along the beam pipe. Application of these 
cuts yields 145,000 t~ — > 7r~7r~7r + j/ T events. 
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Figure 1. E miss /E cm vs. V m i sa /E cm for ex- 
perimental data (left) and Monte Carlo r-events 
(right), where E m i SS is the missing energy and 
Pmiss ^ s ^ ne nrissing transverse momentum in an 
event. This 2-dimensional cut is used to reject 
events with missing tracks and/or energy along 
the beam axis. 



The 7r 7r 7r + effective mass spectrum resulting 
from these cuts is shown in Figure |. The back- 
ground, indicated by the shaded histogram, which 
amounts to about 11% of the events in the spec- 
trum. Monte Carlo studies show that this back- 
ground is from other r-decays: ~ \ comes from 
r-events with missing 7r°('s) and ~ i is from r- 
events with > 1 charged K. Monte Carlo studies 
indicate that backgrounds from multihadronic qq 
events, 77 events, and QED events are negligible. 

3. FIT 

Although it is desirable to characterize the 
structure of the ir~ir~ir + system without model- 
dependent assumptions, this is very difficult in 
practice. In a model-independent approach one 
would assume a general form of the hadronic cur- 
rent that allows the separation and analysis of 
the complex magnitudes of the axial vector and 
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Figure 2. tt + tt + tt~ effective mass spectrum. The 
circles are the experimental data and the his- 
togram is the spectrum resulting from fitting the 
data sample as described in the text. The shaded 
histogram is the background from other r decays. 



pseudoscalar components in bins of Q 2 , s±, and 
S2, where Q 2 is the invariant mass-squared of 
the 7r~7r~7r + system, s\ is the invariant mass- 
squared of the 7r + 7r~ combination that yields the 
larger value of mass-squared, and S2 is the invari- 
ant mass-squared of the other it + it~ combination. 
There would be no assumptions regarding reso- 
nances in the decay. However, in certain chiral 
limits, scalar effects can be induced in the cur- 
rent from the a^(1260) and non-resonant contri- 
butions in the current and thus fake pseudoscalar 
effects ||. Secondly, simply binning in three di- 
mensions is problematical: relatively small bins 
are required to study the r decay structure and 
this leads to some bins with low statistics that 
may be crucial in resolving complicated resonance 
structure. By using a model-dependent fit, one 
can make reasonable assumptions about the reso- 
nant structure and study the likely existence of 
each resonance modeled while using the entire 
dataset in the fitting and without the need for 
binning. 

Using an unbinned maximum likelihood fit, the 
substructure in t~ — » 7r~7r~7r + ^ T decays has been 
fit with the set of resonances given in Table [l| For 
comparison purposes, the current Particle Data 
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Table 1 

Masses and widths of the resonances used in the fits to t~ — » ir~Tr~ir + v T and t~ — * n n ir~v T decays. 
For comparison the current Particle Data Group H masses and widths are listed. 



Particle 




M (GeV) 


T (GeV) 


M PDG (GeV) 


Y PDG (GeV) 


oi(1260) 


1-(1++) 


1.230 


0.400 


1.23 ±0.04 


0.25-0.60 


tt(1300) = tt' 


i-(o-+) 


1.300 


0.300 


1.3 ±0.1 


0.2-0.6 


p(770) 


!+(!-) 


0.774 


0.149 


0.7711 ±0.0009 


0.1492 ±0.0007 


p(1450) 


!+(!-) 


1.370 


0.386 


1.465 ±0.025 


0.31 ±0.06 


/ 2 (127Q) 


0+(2++) 


1.275 


0.185 


1.2754 ±0.0012 


u.iooi_ 0026 


/o(600) = a 


0+(0++) 


0.860 


0.880 


0.4-1.2 


0.6-1.0 


/o(1370) 


0+(0++) 


1.186 


0.350 


1.2-1.5 


0.2-0.5 



Group |6| masses and widths are listed. For the 
Jo (600) and /o(1370) resonance parameters given 
by the Tornqvist unitarized quark model have 
been used. The same resonances and parameters 
have been used in the analysis of r~ — > 7r 7r 7r~^ r 
events [|| selected from the same CLEO II set of 
e + e~ — > t + t~ events. 



ture previously reported for the CLEO II t~ — > 
7T 7r°7r _ t/ r || data, which are based on ~ 15,000 
lepton-tagged events. The fit results on these 
events are given in Table [j| for the reader's conve- 
nience. Both analyses find that axial vector am- 
plitudes with isoscalars, especially a7/(1260) — > 
/o(600)7r - , are large. 



4. RESULTS 

The results of the fit are given in Table | 
where the fti are complex coupling constants. The 
7r~7r~7r + invariant mass spectrum resulting from 
the fit is shown as the histogram in Figure ||. 

As expected, the 0^(1260) — > p°n~ S'-wave de- 
cay dominates. There is a significant amount of 
p-wave ai(1260) decay into isoscalars [/o(600), 
Jo (1370)] plus a pion. r-decay via the pseu- 
doscalar 7r~(1300) is small, although its inclu- 
sion improves the fit. Dalitz plots for sequential 
7r - 7r~7r + invariant mass regions are shown in Fig- 
ure U for data and the fit. 

The isoscalar sub-resonances in the decay of the 
ar(1260) [/ (600), / 2 (1270), and / (1370)] play 
a large role in the fits. To demonstrate the need 
for these isoscalars, the data were fitted without 
them. Projections of the Dalitz plat without the 
presence of these isoscalars are shown in Figures || 
and J7[ For the projection onto the si axis (Fig- 
ure g), the lack of all three isoscalars results in 
a poor fit for Q 2 > 2.25 GeV 2 . For the projec- 
tion onto the S2 axis (Figure 0), the presence of 
/o(600) is needed for all regions of Q 2 , especially 
those of lower Q 2 values. 

These results support the resonance substruc- 
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Table 2 

Preliminary results based on ~ 145,000 t~ — > 7r~'7r _ 7r + i' r events. The uncertainties are statistical only. 
The branching fractions do not sum to 100% due to interferences among the amplitudes. The phases are 
relative to the dominant ai(1260) — > pir s-wave decay. 



Branching fraction (%) lZe(f3) Xm(f3) 



Cli — > piT 


s- 


wave 


69.77 


1.00 


0.0 


ai -> p(1450)7r 


S- 


wave 


1.58 ±0.25 


0.05 ±0.02 


-0.25 ±0.02 


ai — ► (07T 


D 


-wave 


1.68 ±0.16 


0.61 ±0.03 


0.35 ±0.05 


ai ->■ p(1450)7r 


D 


-wave 


2.77 ±0.29 


-1.60 ±0.11 


-1.21 ±0.10 


ai - / 2 (1270)7r 


P- 


-wave 


0.68 ±0.08 


-0.12 ±0.07 


0.96 ±0.06 


ai - /o(600)tt 


P 


-wave 


48.40 ± 1.61 


1.99 ±0.06 


2.70 ±0.06 


ai - /o(1370)7r 


P- 


-wave 


11.51 ±0.82 


-0.04 ±0.04 


-1.15 ±0.04 


7T(1300) /97T 


P- 


-wave 


0.01 ±0.02 


0.001 ±0.002 


-0.002 ±0.001 


tt(1300) -> /9(1450)7r 


P- 


-wave 


0.02 ±0.03 


-0.005 ±0.010 


-0.013 ±0.009 


tt(1300) -> /o(600)7T 


S- 


wave 


0.13 ±0.07 


-0.005 ±0.003 


-0.008 ±0.002 



Table 3 

Results for the substructure fit to ~ 15, 000 lepton-tagged t~ — > 7r°7r 7r~z/ T events from the same CLEO II 
e+e~ — > t+t~ dataset. The branching fractions do not sum to 100% due to interferences among the 
amplitudes. 



Significance Branching fraction (%) 



phase ip/ir 



pir 

p(1450)7r 
pir 

/9(1450)7T 

/ 2 (1270)tt 

/o(600)tt 

/o(1370)tt 



.S-wave 
S- wave 
D-w&ve 
P-wave 
P-wave 
P-wave 
P-wave 



1.4(7 

5.0ct 

3.1(7 

4.2cr 
8.2ct 
5.4a 



68.11 
0.30 ±0.64 
0.36 ±0.17 
0.43 ±0.28 
0.14 ±0.06 
16.18 ±3.85 
4.29 ±2.29 



1.00 
0.12 ±0.09 
0.37 ±0.09 
0.87 ±0.29 
0.71 ±0.16 
2.10 ±0.27 
0.77 ±0.14 



0.0 
0.99 ±0.25 

-0.15 ±0.10 
0.53 ±0.16 
0.56 ±0.10 
0.23 ±0.03 

-0.54 ±0.06 
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Figure 3. Dalitz plots for sequential bins in Q 2 = M^_ 7r _ jr+ . Q 2 is in units of (GeV) 2 . si and S2 are 
the two possible ir + n~ invariant masses squared; si is the larger of the two possible 7r + 7r~ invariant 
mass-squared values. 
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Figure 4. Projection of the Dalitz plot onto the 
si = M^ +7r _ axis for various bins in M| w . The 
points are the experimental data and the his- 
togram is the result of the fit. The shaded regions 
represent the background. 



5. SUMMARY 

This work is a new measurement of substruc- 
ture in t~ — > 7r~ 7r~ n + v T . A model-dependent un- 
binned maximum likelihood fit on 145,000 events 
was performed using an axial vector component 
[ai(1260) primary resonance and several subres- 
onances, including isoscalars] and a pseudoscalar 
component [7r(1300) plus several subresonances] . 
As expected, this decay is dominated by t~ — » 

(1260)^ r with subsequent a\ — > p°Tr~ S-wave 
decay. However, it is also apparent that other 
a 1 "(1260) decay modes are present. In particular, 
0^(1260) decay into isoscalars [/ (600), / 2 (1270), 
and / (1370)] is significant. The aj"(1260) -»• 
/o(600)7r _ is especially prominent. These results 
are consistent with the earlier CLEO II results 
for t~ — > 7r°7r°7r _ ^ r ||. Both analyses find that 
axial- vector amplitudes with isoscalars, especially 
a j" (1260) — » /o(600)7r _ , are large. Regarding 
non-axial-vector contributions, inclusion of the 
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Figure 5. Projection of the Dalitz plot onto the 
s 2 = M^ +7r _ axis for various bins in M| w . The 
points are the experimental data and the his- 
togram is the result of the fit. The shaded regions 
represent the background. 



J p = 0~ 7r~(1300) improves the fit to the data 
in the expected regions, but the need for it is not 
compelling. It is hoped that the future CLEO-c 
dataset will clarify some of these issues. 
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